Theoretical study of ionization of an alkali atom adsorbed on a metal surface by laser 

assisted subfemtosecond pulse. 
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The first numerical simulation of the process of ionization of an atom adsorbed on a metal surface 
by the subfemtosecond pulse is presented. The streaking scheme is considered, when a weak sub- 
femtosecond pulse comes together with a strong IR pulse with a variable delay between them. 
The problem is analyzed with numerical solving the non-stationary Schroedinger equation in the 
cylindrical coordinate. The results obtained are compared with ones in the gas phase. We show that 
the surface influences the DDCS, but the observation of this influence, beside the trivial polarization 
shift of the energy of the initial state, requires a quite high experimental resolution. 
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Development of sub-femtosecond experimental tech- 
nique initiated in 2001 with the paper [l[ , providing a 
tool for monitoring electrons with sub-atomic resolution 
in both space and time, creates principally new oppor- 
tunities for real-time observation and control of atomic- 
scale electron dynamics 0]. With keeping track on the 
temporal evolution of the outgoing electron wavepackets, 
this technique gives direct time-domain insight into var- 
ious aspects of interaction of an electron with very fast 
relaxing many-electron systems. The attosecond streak- 
ing spectroscopy has already enabled observation of the 
decay of an inner-shell vacancy through Auger relaxation 
in isolated atoms in the gas phase 0] and lead to creation 
of attosecond chronoscopy Extension of the subfem- 
tosecond technique to solids is of paramount importance. 
This is because a vast number of transient electron phe- 
nomena (e.g., charge screening, collective dynamics in 
metals and semiconductors, charge transfer in hostguest 
systems and other processes linked with electronelectron 
interaction) at surfaces and in solids does evolve on a few 
femtosecond - subfemtosecond time scale. 

Although a huge number of papers on attophysics 
[j| are already published, a very few proof-of-principle 
experiments have been performed as yet. The time- 
dynamics has been observed with an ingenious 'streaking 
camera' @, when the IR pulse producing the atto-pulse 
by the high-harmonics generation, is simultaneously used 
for transposing the electron ejection time into the elec- 
tron energy spectrum. In parallel with the experimental 
studies, development of theory, which should provide full 
understanding of the experimental findings, is of extreme 
importance. In general, the time-independent theoretical 
paradigm must be supplemented by proliferating time- 
dependent theories directly describing attosecond exper- 
iments in the time-domain. At the current level of devel- 
opment of computational technique it becomes possible 



to build up rather realistic models of the phenomena with 
atoms in the gas phase 0] which are feasible for full, step 
by step, numerical control of the results. Meanwhile, 
comprehensive understanding of dynamic behaviour of 
electrons in condensed-matter systems is much more im- 
portant for promoting experimental investigations and 
their application in nanotechnology. The first experimen- 
tal extension of the attosecond streaking technique to the 
processes with the solid surfaces has been recently per- 
formed 0. This is a very important proof-of-principle 
experiment which has shown that it is possible to ob- 
tain direct time-domain access to charge dynamics at the 
near-surface region of metals by probing photoclectron 
emission from a single-crystal. Theoretical simulations 
of this experiment have been presented in 

Further development should be directed to investiga- 
tion of various processes with surface impurities, includ- 
ing charge transfer, charge screening, image charge cre- 
ation and decay, etc. The simplest object for such stud- 
ies is a single adsorbate on solid surface. In this case, 
ejection of an electron from an internal shell of the ad- 
sorbate produces abrupt change of charge distribution 
external for the solid. Observation of the energy and 
angular distribution of the ejected electrons can provide 
information regarding the relaxation of the electron sys- 
tem at the surface of the solid. Current state of the 
computational physics possesses a potentiality for pro- 
viding some data on the possible effects, which can be 
of importance in planning of the experiments. Here we 
present the first theoretical simulation of streaking exper- 
iment with an atom adsorbed on a metal surface, restrict- 
ing ourselves with consideration of the problem within 
the static one-electron approximation and do not touch 
description of the relaxation processes. We consider a 
solitary Na adsorbate which is placed at Z a d = 2.5 a. u. 
from the image plane of Ag(100), which is at the dis- 
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tance Zi m = 2.064 a. u. from the outermost layer of 
the metal toward the vacuum. The adsorbate is initially 
ionised to Na + state. Then it is ionized once more from 
the 2p-shell by the XUV atto-pulse (duration of about 
200 as, we consider the pulses with photon energies 50 
and 90 eV). The ejected electron moves in the field of the 
Na ++ ion U-2(r), in the field of the lattice, U s (z), and in 
the field of a relatively strong IR pulse with ujjb, = 1.6 eV, 
which is damped in the metal very efficiently. The elec- 
tron ejected into the metal has a rather short mean free 
path (MFP) in the metal (A/ r w 10 a. u.). Both the light 
pulses are assumed to be polarized perpendicular to the 
surface. In this case, the component of the angular mo- 
mentum of the ejected electron perpendicular to the sur- 
face is a good quantum number. The ionization can pro- 
ceed from both a and 7r states, strongly localized in the 
vicinity of Na nucleus. Studying the electron ejection 
close to the normal to the surface, we presently restrict 
ourselves with consideration of only the a-case. 

The basic Schroedinger equation governing the time 
evolution of the system reads: 

*§t*(f,t) = (H(r - R at ) - i 1 -^ - E)W,t) + 

(U surf {z)+E IR {t)z)^{r,t) + (1) 

^Exuv(t)z^ 2p (r- Rat) 

The Hamiltonian H(r — R a t) describes the interaction 
of the electron with the Na ++ core, the correspond- 
ing potential has been computed with the Hartry-Slater 
approximation. The (pseudo)potential U s (z) in Eq.(l) 
mocks the interaction of an electron with the lattice in a 
finite metal. For this potential, we use the parametriza- 
tion by Chulkov et al [T3| with the parameters for 
Ag(100) surface. The damping function 7(2:) in Eq.(l) 
takes into account ineleastic electron-electron collisions; 
it is non-zero at z < 0, where 7(2) for the electron with 
energy E in the bulk is approximated as y/2E/Xf, with 
Xf being electron clastic MFP near the metal surfaces. 
It weakly depends on the electron energy, being close to 
5 A at a few tens eV range We have set A/ = 10 a.u. 

Since the XUV field is rather weak, the interaction 
between the active electron and XUV atto-pulse is de- 
scribed within the rotating wave approximation (RWA, 
for some details, see Q). The quantity E in (1) is a 
sum of the electron energy in the initial state, — |e2p|, 
and the carrier frequency of the XUV pulse, lux ■ The 
wave function of the active electron in the initial state, 
^2 P {t — Rat), has been precalculated to take into account 
the polarization of this state by the surface potential. 
The wave function is slightly changed, but this interac- 
tion strongly influences the energy of the initial state: 
£2 P = —1.3922 a.u., while the energy of corresponding 
state for a single ion Na + within our model is £i v = 
— 1.6895 a.u. (the experimental value of this energy is 



£ 2p P = —1-738 a.u..) The shift of the energy of the adsor- 
bate state is close to the classical image shift +3 /AZ a d = 
0.3 a.u., for an electron localized on the adsorbate center. 
The source term l/2E X uv(t)z^> 2p (r — R at ) in Eq.(l) de- 
scribes ionization of the initial electron state by the XUV 
pulse with the envelope Exuv{t). The interaction with 
the electromagnetic field is taken in the length form, the 
factor 1/2 is due to the RWA used. For the envelope, we 
have used the form ex(t) = exp(— (i — t delay) 2 /t x ) with 
fx = 125 as, FWHM = t x = 0.21 fs. The E L (z) in 
Eq.(l) describes the NIR laser electric field: 
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with the screening length £ = A a.u. [12]. The param- 
eter Zi m in Eq.(4) is a position of the image plane. It 
enters the parametrization of the pseudopotential U s (z) 
(see 0, Eq.(2)-(5)), for Ag(100) z im = 2.064 a.u.. As in 
the experiment [8j , we assume that the streaking field in- 
cidents on the surface at the Brewster's angle (8b ~ 77.7° 
for silver). Thus, the normal to the surface component of 
the IR field in the metal is about 20 times weaker than the 
incident field. For the envelope of the NIR pulse we use 
the form e L (t) = 0.5£ {1 - cos[7rf/T L ]}, < t < 2t l 
with the FWHM t l = 5 fa and lu l = 1.6 eV, = 
in Eq.(l). The field strength £q is conventionally re- 
lated to the intensity of the XUV pulse W , we set 
W = 10 12 W/cm 2 

At variance with the case of a single atom in the gas 
phase, when computation can be performed with the 
spherical coordinate system Q, in the present case we 
are to take into account the strong field by the nucleus 
in the cylindrical coordinate system that causes substan- 
tial technical problems. Briefly, the nonuniform meshes 
are used in both the p and z variables, both the meshes 



are almost cubic (Az 



) at the center of the adsor- 



bate atom (see [l3j]). This allows us to take into account 
the strong core Coulomb potential with high accuracy. 
The computations have been performed with the split- 
propagation scheme. Strong concentration of the mesh 
points on the adsorbate core requires some changes in the 
Crank-Nicolson algorithm (the computational details will 
be presented elsewhere) . The mesh has covered uniformly 
the cylinder z g [-80, 376], p € [0,450] . The time step 
equals to 0.03 a.u. The outgoing wave asymptotic con- 
dition has been provided by the artificial adsorbing po- 
tentials at the edges of the mesh. Propagation has been 
performed till tfi n = txuv + 200 a.u.. Taking into ac- 
count that the electron velocity in vacuum is in 1-2 a.u. 
range, during this time the essential part of the wave 
packet passes through the registering screen (see below) , 
which has been placed at zq = 290a. u. Another problem 
in computational scheme with the cylindrical coordinate 
system is related to a procedure of extraction of observ- 
able quantities from the results of our time-dependent 
computations. While in the case of a solitary atom one 
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can completely describe the spectrum of the eigenstatcs 
of the ion in a free-field case 0], presently construction 
of such states is not feasible. Here we put forward a dif- 
ferent method develops au the method described in 14 1 . 

In a far-zone, where only a time-dependent uniform 
electric field polarized along the z— axis acts on the elec- 
tron, the basis functions, normalized at large time over 
the momentum scale, can be written as: 

^ (p, z, t) = cxp [t(Q(t) + P(t)z + qp)] . (3) 

The Schroedinger equation 

d r 1 i 

l gl^,k„ (Pi z i t ) = [ - g A + £, (*) Z J V'g.fc, (p, 2, t) 

leads to the explicit representation for Q(t), P(t): 

/oo 
cItE (t); 

Q(t) = -(k 2 z + q 2 )t/2 + k z Z x (t) + Z 2 (t)/2; (4) 

/OO pOO 
dTA(r); Z 2 {t) = J dr A 2 (t) 

The basis functions (3) read 

Vfr,fc,(£*M) = ( 27r )3/2 6X P ( l -(^ + 3 2 )*/ 2 

+fc,(« + Zi(t))+A(t)z + Z a (t)/2 + g?]). (5) 

These functions merge with the conventional plane waves 
at t — +oo. The expansion of any wave function 

/>oc 

&(p,Z,t)= dk z d 2 q f(q,k z )ip^ kz (p,z,t), (6) 
Jo 

which contains only outgoing waves, can be rewritten as 
*(p,z - Zi{t),t) = e 4Mt)(^z 1 (t))+z 2 (t)/2] 



dk z d 2 q 



e*[-«+« a )*/3+*.*o+^ mkz) . ( 7 ) 

'o (27r) J /^ 

Thus, the amplitude f(q,k z ) is the Fourier transform of 
the wave packet value at the moving screen z(t) = zq — 
Z 1 (t), where z = z(oo), e = {k 2 z + q 2 )/2 : 

poo 

f(k z ,q)^(27r)-^ 2 k z e^ Z0 I dte ^t-A(t)z(t)-z 2 (t)/2] 

— OO 

d 2 pe-^(p,z(t),t) (8) 
Finally, the double differential cross section reads 
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Being compared with the method of extraction of in- 
formation from the time-dependent computations used 
before , the present approach has a noticeable and uni- 
versal advantage. Within the previous approach, we are 
to perform the computations till the very IR pulse termi- 
nation and, at the same time, are to guarantee that the 
essential part of the wave packet is still in the mesh. This 
requires using a very large mesh in a long IR pulse case 
and, while using the length gauge, this causes at the same 
time decrease of the time step. Thus, within the method 
used previously, the computations for IR pulses longer 
than 10 fs may become unfeasible. Within the present 
approach, one should keep track on the wave packet prop- 
agation only till it passes the 'screen'. The only restric- 
tion on the mesh size is that the moving 'screen', while 
its oscillations, should not enter the region where the ab- 
sorbing potential is operative. Also, the screen should 
all the time be in the region where all fields, beside the 
uniform electric field, can be neglected. This allows one 
to deal with rather long IR field. However, we note that 
the current approach allows to compute the DDCS for 
the electron ejection in a restricted sector in forward and 
backward directions. 

Some results of our computations are given in Fig.l. 
These computations are performed for a fixed delay be- 
tween the center of the XUV atto-pulse and the upset 
of the IR pulse. For such a delay, the IR field at the 
center of the XUV pulse is rather small and therefore 
the vector-potential reaches its extremum, Aq = —0.128 
a.u. Within the conventional theory of the streaking ef- 
fect, the energy of the DDCS maximum is equal approx- 
imately to e max = e 2p + uj + v /2(e 2p + u)A + A\. The 
comparison of the maxima positions is in the table 1 . The 
correspondence is very good except the case of low fre- 
quency uj = 45 eV, when the whole structure of DDCS is 
strongly perturbed with the surface (compare Figs. la and 
lg). Although this frequency range is out of the current 
experimental studies [2| , it could be of interest because at 
such low frequencies the ejected electrons have the veloc- 
ity comparable with the Fermi velocity in the substrate 
and therefore in this case the relaxation of the electron 
system in the metal after the ejection electron from the 
internal shell of the adsorbate could reveal itself most 
noticeably. Clearly is seen an effect of the conventional 
polarization shift in the energy of the initial state due 
to the interaction of the electron with the image charges 
induced in the metal. Currently the experimental studies 
with atto-pulses are concentrated in the frequency range 
about 90 eV. We have obtained a principally noticeable 
effect of the surface on the DDCE in this frequency range, 
but the observation of this influence requires rather high 
experimental resolution. 

In conclusion, a method for numerical simulation of 
DDCS of electron ejection from an internal shell of an 
atom adsorbed on a metal surface by attosecond XUV 
pulse accompanied with strong IR laser field is developed. 
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FIG. 1: The results of computations of DDCS So jdedQ, as 
functions of the final electron energy (e) and angle between 
the final electron momentum and the normal of the surface 
(9) for different frequencies of the XUV pulse (ujxuv). In 
Figs.a-f the DDCS for ionization of the adsorbate are given, 
in Figs. g,h the results for ionization of a free Na + ion. 
a: ujxuv = 45 eV; b: ujxuv = 75 eV; c: ujxuv = 80 eV; 
d: uixuv = 85 eV; e: ujxuv = 90 eV; f: ujxuv = 95 eV; g: 
ujxuv = 53.1 eV; h: ujxuv = 98.1 eV 



TABLE I: Comparison of the energies (in eV) of maxima 
DDSC obtained in computations e ma x and e e // obtained 
within the conventional theory of streaking. 
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It is shown that the surface produces a noticeable effect 
in DDCS. The method can be applied for numerical sim- 
ulation of the effects of relaxation in the electron system 
of substrate caused by the abrupt change of the external 
charge distribution due to absorption of an attosecond 
pulse. 
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